Dynamics of disordered vortex matter in type II superconductors. 
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Dynamics of homogeneous moving vortex matter is considered beyond the linear response. The 
framework is the time dependent Ginzburg - Landau equation within the lowest Landau level ap- 
proximation. Both disorder and thermal fluctuations are included using the Martin-Siggia-Rose for- 
malism . We determine the critical current as function of magnetic field and temperature J c (B,T). 
The surface in the J — B — T space defined by the function separates between the dissipative moving 
vortex matter regime (flux flow) and an amorphous vortex "glass". Both the thermal depinning and 
the depinning by a driving force are taken into account. The static irreversibility line, determined 
by J c (B, T) — is compared to experiments in layered HTSC and is consistent with the one ob- 
tained using the replica approach. The non-Ohmic I-V curve (in the depinned phase) is obtained 
and compared with experiment in layered superconductors and thin films. 

PACS numbers: 74.40.+k, 74.25.Ha,74.25.Dw 



In type II superconductors the magnetic field pene- 
trates the sample in a form of Abrikosov vortices, which 
interact strongly creating an elastic "vortex matter" . Im- 
purities greatly affect the thermodynamic and especially 
dynamic properties of the vortex matter. In high T c su- 
perconductors (HTSC) thermal fluctuations also influ- 
ence the vortex matter, either directly by melting the 
vortex lattice into a liquid or by reducing the efficiency 
of the disorder (depinning). As a result of the delicate 
interplay between disorder, interactions and thermal fluc- 
tuations even the static H — T phase diagram of HTSC 
is very complex and is still far from being reliably de- 
termined. Once electric current J is injected into the 
sample, the dynamical phase diagram should be drawn 
in the three dimensional T — H — J space, which makes 
the analysis even more complicated. Generally there are 
two phases, the pinned phase in which the linear resis- 
tivity vanishes and a phase in which vortices can move 
due to Lorentz force and thus a finite resistivity appears. 
The surface separating the two phases is determined by 
the critical current (neglecting very small creep effect) as 
function of magnetic field H and temperature T. The 
intersection of the surface with the H — T plane gives the 
static irreversibility line. 

Theoretically two major simplifications are generally 
made. In London approximation (valid far from Hri) 
the vortex matter behaves as an array of elastic lines [1| . 
An alternative simplification to the vortex matter, valid 
far from H c \, where the magnetic field is nearly homo- 
geneous due to overlaps between fields of the vortices, 
is the lowest Landau level (LLL) approximation 0, 
The original idea of the vortex glass and the continuous 
glass transition (GT) was studied in the static frustrated 
XY model 0, using RG, variational methods and nu- 
merical simulations |fj. It was found that in this model 
the conductivity exhibits the glass scaling. In analogy 
to the theory of spin glasses, in this model the replica 



symmetry is broken when crossing the GT line. More re- 
alistic model, the elastic medium of interacting line-like 
objects subject to both the pinning potential and the 
thermal bath, was treated using the Gaussian approxi- 
mation 0, and RG [f| . The dynamics in the presence 
of thermal fluctuations, within the latter model, can be 
simulated using the thermal bath Langevin force . 

The irreversibility line (along with other properties) in 
2D and 3D disordered vortex matter was found by ap- 
plying the replica method to the Ginzburg-Landau (GL) 
model [lfj, (supersymmetry was also employed to de- 
scribe the effect of columnar defects in layered HTSC 
[13). The irreversibility line of YBCO and a 2D or- 
ganic superconductor were in good agreement with ex- 
periment. Dynamics in the presence of thermal fluctua- 
tions and disorder is generally described using the time 
dependent Ginzburg - Landau (TDGL) model, in which 
the coefficients have random components . This model 
was studied by Dorsey, Fisher and Huang in the 
homogeneous (liquid) phase using the dynamic Martin- 
Siggia-Rose (MSR) approach 0] . They obtained the ir- 
reversibility line and claimed that it is inconsistent with 
experiments in YBCO. 

It is the purpose of this paper to study the dynamics 
of vortex matter within the TDGL model beyond linear 
response using the dynamical approach of We cal- 
culate the I-V curve in the homogeneous flux flow phase. 
The critical surface in the three dimensional T — H — J 
space, separating the pinned and unpinned phases is ob- 
tained. The static GT line (zero current) coincides with 
the one obtained using the replica method A re- 

lation between the dynamical and the replica methods, 
which in our mind is crucial for understanding the na- 
ture of any glass transition, is discussed. Comparison 
of the irreversibility line, critical current and resistivity 
with experimental results in layered superconductors and 
thin films is made. 
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Our starting point is the TDGL equation in the 
presence of thermal fluctuations on the mesoscopic scale 
represented by a white noise £: 



fc2 

n 7 
4m* 



dtp 



(1) 



where m* is the effective mass of Cooper pair and 7 is 
the inverse diffusion constant. The covariant time deriva- 
tive is D T = -Jpp + ^|-$, where $ is the scalar poten- 
tial describing the driving electric force. The variance 
of the thermal noise £ determines the temperature via 
the fluctuation-dissipation relation: (£ (x, t) (* (y, r')) = 

6 (t — t') 5 (x — y) $^*T. The static GL free energy in- 
cluding the AT C disorder is: 



F= d 



2m 



W| 2 -a'(l + [/)h/f + -|V| 4 , (2) 



where the disorder variance is expressed via dimension- 
less pinning strength n and the coherence length £: 
(U (x) U (y)) = 5 (x — y) £ 2 n. The covariant derivative 
D = V + j— A describes the magnetic field, and the co- 
efficients in Eq.||2J) are related to the coherence length 
and magnetic penetration depth A, namely a' (T = 0) = 



and b' 



written as H tb 



2-Kh 2 \ 2 e* 2 
£ 2 c 2 m* 2 

= a'Uip 



The TDGL equation can be 
b'\ip\ 2 ip + (, where the (non 

I 



Hermitian) linear operator is defined by H = j^D T — 



-tD 2 



Several assumptions (identical to those used 



in 16] and major parts of |3j) are made. In strongly type 
II superconductors, where k >> 1, magnetic and electric 
fields are very homogeneous, since fields of vortices over- 
lap. Therefore the Maxwell equations for electromagnetic 
field are not considered. The Landau gauge with vector 
potential A = (— By, 0,0) and scalar potential $ = Ey 
is used. Temperature, current and magnetic field should 
be close enough to the dynamical phase transition line 
H C 2(T, J) in order to apply the GL approach. 

In order to perform the averaging over both the ther- 
mal fluctuations and disorder in the dynamical situation, 
we adapt the MSR formalism to the present case. 
The dynamical "partition function" is defined as a func- 
tional integral over the order parameter ip and an addi- 
tional " ghost" field cf> (originating from integral represen- 
tation of delta function): 



Z = J DtfjD(f>exp {-Amsr 



[im, in)- 



(3) 



Although the formalism can be applied to both the 3D 
and the 2D cases, to simplify the discussion, wc consider 
a 2D superconductor of thickness L z . The MSR dimen- 
sionless "action" A, with disorder averaged out is: 



i-MSR 



'' T 



(j>* t Hip rt + cc 



: ?i 2 7 
'2m* 



t>r,t\ 2 + b'\lpr,t\ 2 y>*,ti>r,t + Cc] + irb' J g [</>* t Vv,t + cc] [0* s VV,« + Cc] 



(1-t) 2 



Gi ^ \ (H^j 



where r — — „ , — — : — mm _ - • r— -. 

2ir 2 y r 2Gi t 2 ^ 7rc 2 h 2 Lz 

are the disorder dimensionless parameter and the 2D 
Ginzburg number respectively. In high enough magnetic 
field, (j) and if) can be expanded in a LLL basis (right 
eigenfunctions of the operator H with lowest eigenvalue) 
<Pku — exp[z (u)t + kx)] exp — | (?//£ — fc£/6 + ro/6) 2 . 

The dimensionless magnetic field and velocity are b — 
B/B C 2 and v = e*jE£ i 3 /4frb respectively ( 7 depends on 
temperature 7 = 70 (1 — t) 11 , for example in BCS r\ = 1, 
70 = 487TK 2 <7 r i/c 2 , where a n is the normal state conduc- 
tivity |l|). The LLL basis will be used from now on. 
The model is still a highly nontrivial field theory and we 
use the Gaussian effective action approach to treat it. 
The approximation was applied to the LLL Ginzburg - 
Landau model in the absence of electric field (using di- 
agram resummation) in Ref. |l3| . leading to identical 
gap equation. Two point Green functions are the cor- 
relator C (r,t,r' ,t') = lij}*, t /Vv,t)j the response function 
R (r, t, r' , t') = t/"0r,t) and the auxiliary field cor- 



relator B (r, t, r' , t') — (<t)Z, t ,(j)r.t)- In a homogeneous dy- 
namical phase (stationary flow) |l7j , the correlators time 
dependence is just on the time difference t — t' . Within 
LLL we find a solution, which does not depend on mo- 
mentum k : C {k,t,k' ,t r ) = C[t — t'). The variational 
parameters therefore will be R u , and B u , which we 

write in a matrix form q = ( . ^ ] . 

' \ -iR B J 

The Gaussian effective action is: 
A (X J [Trlogg- 1 +Tr(g^g)]+ L J C X (R U -R*J 

UJ.X 

+ -y J [^BujC u - i?.* 2 - i? 2 ] , 

where = 4ntbV2Gi, a h = - (l - t - b - v 2 ) /2 and 

_! = L z £h 2 e v2 l b f -^7^/2 - 2ia h \ 

90 ~ 2 7 r 1 /2 m *T6 3 / 2 I ^ 2 /2 - 2*a h 7 ^ 2 J ' 



FIG. 1: In (a) the imaginary and real parts of the response 
function are shown for arbitrary values of the parameters. 
The correlation function has a peak at u — 0. In (b) it is 
shown qualitatively how the correlation function at zero fre- 
quency diverges as the parameters approach the critical val- 
ues. 



The corresponding " gap equations" are 

(2a h + 6C(t = 0)+ R u - rfli? 2 = 1 (4) 



r lew .o _ n 

1 — r6 \FL,\ 



(5) 



For v = the equations are consistent with those of Ref. 
IfJ. For cj = 0, one obtains 




FIG. 2: The Irreversibility surface in the dimensionless j, t, b, 
parameters space. 



given by r max = which diverges as one approaches 
the critical surface determined next. 

There are two dynamical phases of the system in the 
tree dimensional external parameter space (T, H, J). It is 
more convenient to use instead dimensionless scaled pa- 
rameters (t, b, v) . The critical surface is defined as a set 
of values of the parameters for which the correlator C w at 
ijj = diverges. We will argue later that below this sur- 
face the superconductor acquires certain "glassy" prop- 
erties. According to Eq. JSJ), the correlator at zero fre- 
quency diverges when the denominator vanishes, namely 



-a h + y/aj + Ojl-r) 
9(l-r) 



(6) 



Eqs. (|5I6|) are consistent with the static (t = and v = 0) 
correlator calculated using the replica method |h] . The 
solution for arbitrary lo is: 



Ru) — 
where 



4r9 



_ a h (1 - 2r) + y/gi + (1 - r) 
(1-r) 

The correlator and the response functions versus fre- 
quency are plotted in Fig. ^ The correlator decreases as 
1/ui 2 at large frequencies. The structure of the singular- 
ities in the complex oj plane is as follows. In addition to 
cuts there is a segment of singularities along the positive 
imaginary axis between uj™* = 2~f~ 1 £~ 2 (a v ± V 4r9j . 

This range of frequencies corresponds to a range of re- 
laxation times. Asymptotically the long time dependence 

of the correlator is exponential, C (t — t') cx exp • 

The dominant time scale for relaxation is the longest one, 



l-r9\Rl =0 \ 2 =0. 



(7) 



Using the response function Eq. ©, one obtains the 
critical surface 



4 ( w ) =4r 1 / 2 -2r- 1 / 2 , 



(8) 



where the LLL scaled temperature is ar (v) = 4ai l 9~ 1 < >2 . 
In terms of the original dimensionless parameters t, b, v 
it takes the form: 



\-t-b-v 2 = 9 1 ' 2 r- 1 ' 2 -2r 1 ' 



(9) 



The static irreversibility line (v — 0), which is the in- 
tersection of the dynamical glass transition surface (Fig. 
EJ with the H — T plane, has been observed in great 
variety of type II superconductors in magnetic fields. 
Within the range of applicability of the LLL approxima- 
tion the static irreversibility line in both 2D and 3D was 
already obtained using both the replica method and 
the dynamical approach |l3l| . The original impression at 
\\?\ was that the line is inconsistent with experiments in 
thin films of YBCO, however recently the irreversibility 
line in organic superconductor k type BEDT — TTF, 
was found to be in good agreement with the experiment 
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FIG. 3: The glass line in the b-t plane. The pluses correspond 
to experimental data of |T3| . and the solid curve corresponds 
to the analytical formula Eq. Q, where v = 0. 



FIG. 4: Resistivity as function of the magnetic field. The 
solid curves show the analytical results (Eq. H12I 1. while the 
pluses correspond to experimental results of |l9l| . 



[HI. In Fig. we compare the irreversibility line of 
BSCCO in high magnetic fields []J] with Eq. © for 
J = v = and the following values of the material param- 
eters: H c2 = 195T,T C = 93K, Gi = 0.00044, n = 0.005. 
One would like to parametrize the surface in terms of the 
supercurrent density 0] rather than in terms of the flux 
velocity v. It provides the nonlinear I-V curve: 

2e*vTb 
Js — , tt Ruj=a 

Tic 2 v 



1/2 - Q r (v) + \A| (v) + 16(1- r) 

*\2" isrt J V ") 



16 (1 - r) 



Note that Eq. 1)10(1 is valid only for Js > J c - The critical 
current according to Eq. is: 



Jr. = 



hc 2 Q 1 ' 2 
4e*A 2 £r 1 /2 



l-t-b + 20 1 ' 3 hr 1 ' 2 - r- 1 ' 2 



1/2 



_ (11) 

The critical current of the MoGe films |2fJ as function of 
temperature compares qualitatively well with Eq. 1|11|) 
in the region beyond the peak effect (namely in the ho- 
mogeneous phase). 

In the linear response limit Eq. (|10|l determines the 
conductivity due to Cooper pairs: 



37r3/2 ^V/4 f 1 a ? (°) ~ V q t (°) + 16 (1 - r) 

(12) 

This expression, valid for non-zero electric field, has a 
finite limit when one approaches the critical surface. 
Within the vortex glass theory Q the conductivity di- 
verges, and in 13j it was argued that higher order cor- 
rections to the Gaussian approximation lead to this diver- 
gence. We were unable to confirm that and believe that 
the exponentially small creep may appear when instanton 
effects are taken into account and higher Landau levels 
are added. Comparison of the resistivity in BSCCO 0] 
with the one obtained from Eq. (|12fl (taking into account 



the normal part of the conductivity) is shown in Fig. 0] 
We conclude by discussion of the transition to the glassy 
state. When the distance from the critical surface in the 
parameter space A = ax (v) — a g T (v) approaches zero, 
certain physical quantities diverge powerwise. For exam- 
ple, the relaxation time diverges as r max oc [clt — ciy) 
and C (ui = 0) oc (ay — dj,) 1 ■ Note however that the 
static correlator C (t = t') (proportional to the magneti- 
zation within LLL [lj|) does not diverge at the critical 
surface. Thus it can not serve as an order parameter for 
the "glass" transition, even in the static limit. The com- 
mon wisdom is that replica symmetry is broken in the 
glass (either via steps or via hierarchical continuous pro- 
cess) as in most of the spin glasses theories |2l|. The 
replica method applied to the static LLL model with 
AT C disorder within Gaussian approximation 11] indi- 
cates that there is no replica symmetry breaking in the 
homogeneous phase. However the Edwards -Anderson 
parameter vanish above the GT, while is nonzero below 
it. This is in agreement with the original approach to 
the glass transition of EA (see jl^] for a discussion). The 
results obtained here demonstrate criticality in this case. 

To summarize, using the dynamical approach we ob- 
tained the dynamical critical surface separating the liquid 
and glass phases. In particular the static irreversibility 
line was obtained, and shown to be in a good agreement 
with experiment in BSCCO. The resistivity is also found 
to be in a good agreement with experimental results. I-V 
curve and critical current are calculated beyond the lin- 
ear response limit using the dynamical approach at finite 
electric field. 
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